if (FALSE) {
library(EMCluster, quietly = TRUE)
x2 <- da2$da
emobj <- list(pi = da2$pi, Mu = da2$Mu, LTSigma = da2$LTSigma)
eobj <- e.step(x2, emobj = emobj)
emobj <- m.step(x2, emobj = eobj)
summary(emobj)
ret <- starts.via.svd(x2, nclass = 10, method = "kmeans")
summary(ret)
}
Run the code above in your browser using DataLab